A note on generalized electrodynamics 



S. I. Kruglov 

University of Toronto at Scarborough, 
Physical and Environmental Sciences Department, 
1265 Military Trail, Toronto, Ontario, Canada MIC IA4 



Abstract 

The generalized Maxwell equations with arbitrary gauge parameter 
are considered in the 11 x 11-matrix form. The gauge invariance of such 
a model is broken due to the presence of a scalar field. The canonical 
and symmetrical Belinfante energy-momentum tensors are found. The 
dilatation current is obtained and we demonstrate that the theory pos- 
sesses the dilatation symmetry. The matrix Schrodinger form of equa- 
tions is derived. The non-minimal interaction in curved space-time is 
introduced and equations are considered in Friedmann— Robertson — 
Walker background. We obtain some solutions of equations for the 
vector field. 



1 Introduction 

We have investigated the first-order formulation of the generalized Maxwell 
equations which describe massless vector fields with an additional scalar field 
in [1], [2] (see also |3]). Such a model is not a gauge- invariant and can be 
treated as a Maxwell theory in the definite gauge. Gradient terms were intro- 
duced in Maxwell equations by many authors (see references in [3] and [4]). 
Here, we take into account a gauge parameter which allows us to consider 
different gauges. It should be mentioned that gauge parameter is physical 
value in our scheme, contrarily to classical electrodynamics, which can con- 
tribute to gravity interaction. Therefore, schemes with difference gauges are 
not equivalent each other. As was mentioned in [2j the reason for leaving a 
scalar field in the spectrum is the application of such a non-gauge-invariant 
model in astrophysics. We have stressed [2j that the additional degree, a 
scalar field, can play an important role in the inflation theory of universe. 
Later, authors of the work [6] considered a scalar field of generalized electro- 
dynamics as a source of dark energy. Dark energy is introduced to explain 
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the acceleration of expanded universe at the present time and results in the 
negative pressure. Dark energy interacts only gravitationally representing 
weakly coupled substance. Scenario of inflationary universe [5] allows us 
to understand observable data: our universe is homogeneous and isotropic 
for scales > 100 Mpc (1 pc = 3.26 light years) and expands in accordance 
with the Hubble law. In the chaotic inflation model a massive scalar field 
(quintessence) minimally coupled to gravity is responsible for slow-roll in- 
flation and plays the role of dark energy. But, in this model, the potential 
terms should be fine tuned to have the acceleration of universe at the definite 
time. Another phenomenological way to describe dark energy is to introduce 
a cosmological constant A into the Einstein equation (the term {—Ag^^) in 
the left side of the Einstein equation) . One implies the existence of vacuum 
energy by introducing the cosmological constant. If A > the additional 
cosmological term leads to anti-gravity. But in this case the difficulty arises: 
a vacuum solution is not a Minkowski space-time. In addition, there is no 
physical explanation of a coincidence problem: matter and dark energy den- 
sities possess the same orders of values at the present time and had big 
difference in magnitudes in previous eras. Today, approximately 70% of the 
energy density of the universe is in the form of dark energy, and the rest 30% 
is in the form of non-relativistic matter. Cosmological constant, which gives 
the energy density, remains constant during the expansion of the universe, 
but the energy density of matter and radiation decreases in time. Thus, the 
nature of dark energy is one of the most important problems in astrophysics. 

In the scenario suggested in [6], the time component of a field in gener- 
alized electrodynamics grows in time and becomes dominant explaining the 
acceleration of universe. Therefore, it is of great interest further investiga- 
tion of the generalized electrodynamics where the addition degree, the scalar 
state of the field, can play the role of dark energy. 

The paper is organized as follows. In Sec. 2, the generalized Maxwell 
equations with arbitrary gauge parameter are formulated in the matrix form. 
Matrices of the relativistic wave equation (RWE) obey the generalized Duf- 
fin — Kemmer— Petiau (DKP) algebra. We obtain, in Sec. 3, canonical and 
the Belinfante dilatation currents which are not conserved. The conserved 
modified dilatation current is also found demonstrating the scale invariance 
of the theory of massless fields. We obtain the Schrodinger form of equations 
and the quantum-mechanical Hamiltonian in Sec. 4. The minimal equation 
for the matrix Hamiltonian is found. In Sec. 5, a novel non-minimal inter- 
action in curved space-time is introduced and equations are considered in 
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Friedmann— Robertson— Walker (FRW) background. The solution of equa- 
tions for the time component of the four-potential is found which grows in 
time. We discuss the results obtained in Sec. 6. The quantum-mechanical 
Hamiltonian is found from relativistic wave equation in Appendix A. Start- 
ing with the second-order formulation of the theory, we obtain canonical and 
symmetrical energy-momentum tensors and dilatation currents in Appendix 
B. In Appendix C the quantization of fields is performed in the second-order 
formalism. 

The Euclidean metric is used in Sec. 1-4 and Appendixes A, B and C, and 
four-vectors are = (xm,X4) = {xm,ixo), and xq is a time; Greek letters 
run 1, 2, 3, 4 and Latin letters run 1, 2, 3. We use natural units h = c = 1. 

2 First-order form of equations 

In [7], we considered the general Lagrangian form of massive vector fields. 
For the case of neutral massless vector fields it reduces to 

C = id,A,) (duAp) , (1) 

where 5^i/,o-p = CLSfiuSap + bS^(j6^p + c6^p6a^^. To have the standard kinetic term, 
we put a = —1/2. The Euler-Lagrange equations follow from (1): 

dlA,-2{b + c)d,d^A^ = 0. (2) 

One can see from Eq.(2) that only one parameter (6 + c) remains in the 
equations of motion. It is convenient to choose: c = 1/2, 2b = —C,, where ^ 
defines the gauge. Then Eq.(l) and (2) become 

^ = -\F'^.-lad.A^)\ (3) 

^^.F^u - ^d,d,A, = 0, (4) 

where F^^^, = d^Ay — dyA^ is the field strength. Eq.(4) can be treated as 
the Maxwell equations with the additional gauge parameter. In QED the 
physical values do not depend on the gauge [8], but in our scheme, we expect 
the dependence on ^ because the scalar state presents in the spectrum. In 
P, [2], we have chosen the gauge ^ = 1. At ^ = 0, one arrives at standard 
Maxwell equations. Here we imply that ^ 7^ 0. Introducing notations i^i^y] = 
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{l/K)F^y, ip^ = Af^, ipQ = —{C,/K)d^A^, where k is the mass parameter, 
second order equation (4) can be represented as a system of the first-order 
equations: 

^u^[^,u] + d^ipo = 0, 



(5) 



We note that fields ipA {A = 0, /i, [/iz/]) have the same dimension. Introducing 
wave function = {iIja{x)}, and using elements of the entire matrix 

algebra e^'^ obeying equations: (^'^'^)^^ = ^ac^bd, e^'^e'^'^ = 6bc£^'^, 



Eq.(5) can be written in the first-order matrix form 

'1 



a,yd,y + K l-Ps + Pt 



^!{x) = 0, 



(6) 



where 



a,, 



p. 



_0,0 



p 



1 



(7) 



At the Feynman gauge ^ = 1, Eq-(6) is simplified because Ps + Pt = £^'^ + 
(l/2)e[H-M is the projection operator but + Pt is not. The 11 x 11 

Hermitian matrices obey the generalized Duffin— Kemmer— Petiau algebra 



(8) 



and Pg, Pt are the projection matrices, P^ = P,, P^ = Pt, and extract the 
scalar ant tensor parts of the wave function, respectively. One can verify the 
relations: 



(1) 
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Introducing the Hermitianizing matrix rj pQ: 



the "conjugated" equation reads 



4 4 I 

e ' + e 



m4],[m4] _ [mn],[mn] 

2 



K 



-Ps + Pt 



(9) 



(10) 



where \1/ = \l/+?7 = (^—ipo,ipf^,—ilj[fj.u]j- Matrices and r] satisfy equations: 



can be written as follows: 



In the first-order formalism the Lagrangian 



Otudu + {-Ps + Pt 



^!(x). 



(11) 



Eq.(6) follows from Lagrangian (11) by varying the corresponding action. 
Similar to the Dirac theory, Lagrangian (11) vanishes for fields obeying 
RWE (6). It should be noted that in the second-order formalism, based 
on the equations (4), Lagrangian (3) vanishes only within four-divergence. 
Lagrangian (11), with the help of Eq.(7) becomes 

C = ^od^ip^ - ip^.d^ipo - V'p^mV'M + ^M^mV'p + (^4^1 + \'4^[pij^ ■ (12) 

One may verify that Lagrangian (12) vanishes for fields obeying equations 
of motion (5) and within four-divergence, which does not infiuence on the 
equations of motion, can be represented as 



(13) 



As K has the dimension of the mass, we need to renormalize the fields, and 
under the replacement t/'m iV^/ V^)^fM, Lagrangian (13) coincides with 
(3). Of course, one could define fields ipA according to this replacement from 
the very beginning. It was pointed on the importance of normalization in 
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3 Energy-momentum tensors and dilatation 
currents 

Now, we investigate the scale invariance in the model with arbitrary gauge 
parameter ^. For this purpose, one needs to obtain energy-momentum tensors 
and dilatation currents. The conserved canonical energy-momentum tensor 
(see [2]) is 

^ ^ (14) 

so that dpT^^ = 0. We took here into account that Lagrangian (11) vanishes 
on the solutions of equations of motion. The canonical dilatation current [12] 
becomes 

= x„t;„, (15) 

with its non-zero four- divergence 

The appearance of the gauge parameter ^ here is due to using equations of 
motion (5). The conserved symmetric Belinfante energy-momentum tensor 
is given by [2] 



= 2Ktp[xp]ip[a\] - 2i;^da4'o - 2ipadpi'o 

(17) 



+Sapdp (V'O^/?) - Sapd/B (^1px1p[Xi3]) ■ 



We find non-zero trace of the symmetric Belinfante energy-momentum tensor 
(17): 

T;^^ = 4dp {^0%) . (18) 
A modified dilatation current j2j is given as follows: 

Dj^ = XaT/^a + ^x^[xp] - S^o^Pp. (19) 

One may verify that the divergence of the Belinfante dilatation current 
(19) coincides with the divergence of the canonical dilatation current (16), 
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d^D^ = d^D^^ 7^ 0. As the trace of the Behnfante dilatation current (18) is 
a four-divergence, we may define new conserved dilatation current 

D^ = x„T^^„-4^o^^, (20) 

and d^D^ = 0. Thus, new dilatation current (20) is conserved, and strictly 
speaking, the dilatation symmetry is not broken and the model possesses 
the scale invariance. We have change the conclusion made in [2] about the 
scale invariance because of obtaining new conserved dilatation current (20). 
Relations found in this section are the generalization of formulas obtained in 
[2] on the case of arbitrary gauge ^. 



4 Schrodinger form of equations 

It should be noted that in some cases the Schrodinger equation has ad- 
vantages for the investigation of interacting field problems. To obtain the 
Schrodinger form of equations and quantum mechanical Hamiltonian, we 
have to exclude the non-dynamical components from Eq.(5). For this pur- 
pose, Eq.(5) can be represented as follows: 

^[mn] = - {dmi^n " C^nV^m) , (21) 
K 

Third equation in (21) possesses only spatial derivatives and, therefore, V'[mn] 
are non-dynamical (auxiliary) components. Excluding ip[rnn] from Eq.(21), we 
arrive at the system of equations containing only dynamical components 

1/ ^^^^ 

Let us introduce the 8-component wave function 



$(a;) 



M^) I • (23) 
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Exploring the elements of the matrix algebra, Eq.(22) can be represented as 
follows: 



(24) 



(26) 



Then Eq.(24) takes the Schrodinger form: 

idt^{x)=U^{x), (25) 
where the Hamiltonian is given by 

In the momentum space the Hamiltonian becomes: 

The 8-component wave function (23) describes fields with four spin states 
with positive and negative energies, and there are only dynamical compo- 
nents in the wave function. The matrix Hamiltonian (27) obeys the minimal 
equation as follows: 

[n^ - k^)' [n^ - 2}^) = o. (28) 

There are three eigenvalues of the Hamiltonian: 0, k^, 2k^. The physical 
eigenvalue is k^ so that k\ = \^ . In Appendix A, one can find the Hamiltonian 
(26) expressed through the matrices (7) of equation (6). 
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5 Interaction with gravitation 



5.1 Equations 

Now, we consider the non-minimal interaction of the massless vector field 
with gravity. Some aspects of interactions of vector fields with gravity where 
investigated in [13], [11], [15], [IS], [H] [H], [H], [20], [21]. Let us consider 
the novel action 



S 



-9 



R 



16nG 



(29) 



where G is the gravitational (Newton) constant, A is a coupling constant, and 
are covariant derivatives. In Eq.(29), we have introduced the coupling 
of a scalar curvature with the vector field A^. Within the four- divergence, 
which does not change equations of motion, the non-minimal interaction term 
in action (29) also can be represented as 



int 



OR 

dx^ 



Such term can be added for any vector-tensor theory of gravity but it will 
lead to the higher derivative model. In [22] other couplings to gravity were 
investigated. Let us consider FRW space-time with the fiat spatial part with 
the metric 



gm = 9 = 1, 9u = 922 = 933 = -a{t) , 9 



11 



9 



22 



9 



33 



1 



aity 



(30) 



where a{t) is a scale factor. Nonzero components of the Christoffel symbols 
and curvatures are given by 



T-i _ -pi 

i 10 — J- o 



20 



p3 

30 



a 



11 



22 



pO 
^ 33 



aa 



R 



■00 



-3-, Rik = 6ik (ad + 2d^) , R = -Q 



d 




- + 




a 





(31) 



where H = d{t)/a{t) is the Hubble parameter; R are the Ricci and 

scalar curvatures, respectively. Varying action (29) with respect to vector 
field A^ yields the equations of motion as follows: 



V^F"'' - ^V^V^A^ + XVR = 0. 



(32) 
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For homogeneous electromagnetic fields {diA^ = 0) the fi = component of 
Eq.(32) in FRW background gives 

V^F^" - ^V^V^A" + XV^R = (i° + 3HA° + 3HA^) +XR = 0. (33) 

Taking into account Eq.(31), one obtains from Eq.(33) the equation for Aq: 



Ao + dt(^3HAo-^R'j =0. 



The equation for spatial components are given by 



Am + HAj^ 



0. 



(34) 



(35) 



Variation of action (29) with respect of the metric leads to generalized Ein- 
stein's equation. 



5.2 Solutions 

Let us obtain solutions to novel equation (34). Integrating Eq.(34), one finds 



(36) 



where Ci is the integration constant. Eq.(36) with the help of Eq.(31) be- 
comes the first-order non-homogeneous differential equation: 



Ao + 3-Ao = Ci- — 
a 4 



a 




- + 




a 





(37) 



We can take H = p/t, with p = 1/2 for radiation and p = 2/3 for matter 
eras respectively. In this case, we obtain the solution to Eq.(37): 



U3p- 1) 



(38) 



The last term in Eq.(38) is due to non-minimal interaction introduced in 
Eq.(29). The first term of Ag-component grows with the cosmic time for any 
p [6], and the last term decays. The Aj-component grows but the temporal 
component (Aq) dominates |6]. As was noted in [6j, the term V^A'^ plays 
the role of a cosmological constant during the evolution of the universe. 
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6 Conclusion 



The generalized Maxwell equations with arbitrary gauge parameter are for- 
mulated in the first-order formalism. The gauge U(l)-symmetry of a model 
is broken. As a result, the scalar state of the field presents in the spec- 
trum of the theory. If one introduces the four-current in the right side of 
Eq.(4), then due to the conservation of the current, d^J^ = 0, the equation 
df^diyAi, = holds. It means that the scalar state of the field does not 
interact with charges and currents. But this scalar state can interact with 
gravity via the coupling (29), we have introduced. 

As the matrices of the RWE obey the generalized DKP algebra, one 
can apply covariant methods for finding solutions for definite spin (one and 
zero), spin projections and energy-momentum [I], [2]- Although the canoni- 
cal and Belinfante dilatation currents, found within the first-order formalism, 
are not conserved, we have obtained the conserved modified dilatation cur- 
rent. This demonstrates the scale invariance of the massless fields theory. 
The Schrodinger form of equations obtained possesses some advantages be- 
cause it contains only dynamical components of fields. The found quantum- 
mechanical Hamiltonian can be used for investigation of problems with in- 
teracting fields. RWE, as well as Hamiltonian, are simplified for the choice 
^ = 1 which was used in pQ, [2]. In addition, the commutation relations (C5) 
take the canonical normalized form at ^ = 1. In our opinion the value of the 
parameter ^ = 1 is natural. The consistency of the model will be studied in 
subsequent papers. 

We suggest a novel non-minimal interaction of fields in the FRW back- 
ground. The solution of equations for the time component of the four- 
potential grows in time in the same manner as in [B]. Therefore, the model 
considered has the similar behavior. Although there are difficulties with the 
unbounded Hamiltonian and indefinite metrics (see Appendix C), the model 
has attractive features. We leave the detailed analysis of the model based on 
the action (29) for further investigations. 
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Appendix A 



Let us obtain the Schrodinger equation and quantum-mechanical Hamil- 
tonian from relativistic wave equation (6). One can find from (6) the equation 
as follows: 



ia4dt'^{x) 



aada + K i-Ps + Pt 



(Al) 



Taking into account the relation 04 {a^ — 1) = 0, which follows from algebra 
(8), one can introduce the projection operators: 



A = al = £0,o^^M,A.^^[m4],M]^ n = 1- 



^ P [rnn] , [mn] 



{A2) 



with the properties A = A^, 11^ = 11, All = HA = 0, A + fl = /n (In is unit 
11 X 11-matrix). Operator A extracts 8-dimensional subspace of dynamical 
components and operator 11 acts in 3-dimensional subspace of non-dynamical 
components of the wave function Thus, we introduce dynamical, 4>{x), 
and non-dynamical, x{^)y functions: 

(pix) = A*(x), x{x) = n*(x). (A3) 
After multiplying Eq.(Al) by the matrices and 11, one finds equations 



idt(f){x) — 



Q={al- 1) 



Oiada + « [jPs + Ptj 

Oiada + i^Ps + Pt 



{^{x)+x{x)), 



{(t){x) + x{x)) ■ 



{AA) 



(A5) 



In these equations we imply that the direct sum of functions (f){x) and x(x) 
is ^ {(t>{x) -\- x{x) — ^). One can verify equations 

Qp, + Pt j n = n 0p, + Ptj = n, na„n = o, 

and obtain from Eq.(A5) the function x(x): 



X{.x) = — Uaada(j)ix). 

K 



{A6) 



12 



Excluding the x{^) from Eq.(A4) with the help of (A6) and using the relation 
a^H = 0, we find the Schrodinger equation idt(p{x) = 'H(j){x) , with the 
Hamiltonian 

--a4aallabdadb. {A7) 

K 

Although the matrices are 11 x 11-matrices, the Hamiltonian (A7) acts in 
8-dimension subspace. One can check that Hamiltonian (A7) coincides with 
Hamiltonian (26). Thus, we have obtained here the Hamiltonian in terms of 
matrices of relativistic wave equation (7). It should be noted that Eq.(A6) is 
equivalent to Eq.(21) for non-dynamical components, ipmn- The Schrodinger 
equation with Hamiltonian (A7) does not contain non-dynamical components 
and can be used for solving some problems of interacting fields. 



aada + K + Pt 



Appendix B 

Now, we obtain canonical and symmetrical energy-momentum tensors 
and dilatation currents starting with the second order formulation based on 
the Lagrangian (3). The canonical energy- momentum tensor found from 
equation 

dC 



e 



is given by 



and is conserved: = 0. The trace of the energy-momentum tensor is 

non-zero and reads 



The dilatation current is given as follows |12j : 

dC 



H 
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One can check with the help of equations of motion (4) that the dilatation 
current (B3) is conserved, = 0. Thus, the the scale invariance is valid. 

The canonical energy-momentum tensor (Bl) is not symmetrical. To 
obtain the symmetrical Belinfante tensor, we use the formulas ^12j : 



e 



B 

flQ 



X 



where the matrix elements of the generators of the Lorentz group S^,^, 
Euclidian space-time, are given by 



m 



From Eq.(B4), with the help of Eq.(B5), we obtain 



(55) 



(56) 



so that Xp^y is antisymmetric in indexes and dpd^Xp^y = 0. The sym- 
metric and conserved Belinfante energy- momentum tensor, using Eq.(B4),(B6), 
becomes 



per 



{d^A^f -iApdp{d^A, 



m 



At ^ = Eq.(B7) converts into energy- momentum tensor of electrodynamics. 
Energy-momentum tensor similar to (B7) was found in ^23j by varying action 
on the metric tensor. A modified dilatation current is given by [12] 



D 



B 



r P)^ 4-V 



ill 



where the field-virial V^, in our case, becomes 

Vfj, = Ilf^aAa — Hq,^ (Sq,^)^^ Act = ^^A^ (SqAq) . 

One can obtain the trace of the Belinfante tensor (B7): 
= -^^d, [A, (d^A^)] . 



(58) 



(59) 



(510) 



From Eq.(B8),(B10), we obtain 
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and the modified dilatation current is conserved. Thus, the scale invariance 
takes place with the conserved currents (B3) and (B8). 



Appendix C 

Let us consider the quantization of fields for given Lagrangian (3) (see also 
). Conjugated momenta for generalized "coordinates" A^{x) are given by 

dC dC 

TXo{x) = — T- = -^df^Af,, TTmix) = = + d.^Ao. (CI) 

oAq dAm 

Then the density of the Hamiltonian becomes 

H = TTmA-m - TTqAq - C = A^ (^m + C^m^o) 



One can verify that the equality H = 0^4 holds where the canonical energy- 
momentum tensor 6^^^ is given by (Bl). It should be noted that the classical 
energy S = J Hd^x is not bounded from below and the system is unsta- 
ble. Therefore, we need to introduce indefinite metrics for quantization. 
With the help of standard commutation relations for canonical variables 
[y4^(x, t), ■K^{y, t)] = i5^y6 (x - y), one obtains 



A„(x, t), im(y, t) + 9mAo(y, t) = i6mnS (x - y) 
[Ao(x, t), d^A^iy, t)] = -1<5 (x - y) . 



(C3) 



In the momentum space the real fields read 
1 



k 

where /c^ = — /cq = 0, is the normalization volume. It should be 
noted that the field A^ possesses four independent components: two com- 
ponents are transverse, one component is longitudinal, and one component 
corresponds to the scalar polarization. In [1] , [2] , we found four independent 
solutions to equations of motion in the form of matrix-dyads for the gauge 
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^ = 1. The fields (C4) satisfy commutators (C3) if creation and annihilation 
operators obey the commutation relations as follows: 



am(k),a+(k) 



5rr,.r:.5 ( k - k 



ao(k),a^(k) 



(C5) 



For the Feynman gauge = 1 (which was used in [2]) the RWE (6), (10) (and 
(A7)) are simplified. In this case the "wrong" sign (— ) in the commutator for 
ao(k), cto (k) in (C5) indicates on the necessity of introducing the indefinite 
metrics. For any gauge ^ there are difficulties with the presence of the ghosts 
if one considers four polarizations of the field to be physical. In QED 
the photon fields possess only two polarizations and physical values do not 
depend on the gauge ^ due to the restriction on the physical Hilbert space 

El. 
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